WIND TUNNEL WALL CORRECTIONS FOR ARBITRARY 
PLANFORMS AND WIND TUNNEL CROSS-SECTIONS 



Chester Arthur Heard 



MBnwwuwwt 

|««irOSTOWm»TESCHOOi 




NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

WIND TUNNEL WALL CORRECTIONS FOR ARBITRARY 
PLANFORMS AND WIND TUNNEL CROSS-SECTIONS 

by 

Chester Arthur Heard 
June 1977 

Thesis Advisor: L.V. Schmidt 

Approved for public release; distribution unlimited. 



T178648 




UNCLASSIFIED 



SECURITY classification OF THIS RASE flFhwi Oa(a Snttrtd) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT’S catalog NUMBER 


4. title Cand SuStilia) 

Wind Tunnel Wall Corrections for 
Arbitrary Planforms and Wind Tunnel 
Cross-Sections 


5. TYPE OF REPORT S PERIOD COVERED 

Master's Thesis; 

June 1977 


S. PERFORMING ORG. REPORT NUMBER 


7. AUTMORCa; 

Chester Arthur Heard 


S. CONTRACT OR GRANT NUMBERC*) 


*■ performing organization name and address 
Naval Postgraduate School 
Monterey, California 93940 


10. PROGRAM ELEMENT. PROJECT. TASK 
AREA A WORK UNIT NUMBERS 


II. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


12. REPORT DATE 

June 1977 


13. NUMBER OF PAGES 

75 


U. monitoring agency name & ADDRESS^// diitarwnt irom ControiUng 0iilc9) 


IS. security class, (ol thio riport) 

Unclassified 


I5«. DECLASSI FI cation/ down GRADING 

schedule 


16. DISTRIBUTION STATEMENT (ol thit Rtport) 

Approved for public release? distribution unlimited. 


17. Distribution statement (ot thm mb9trmct 9*%t9rmd Jn 3/ock 20, i/ d//f«r«nf ^rocn R9portJ 


IS. supplementary notes 


19. KEY WORDS ('Continu9 on f9V9f99 9id9 it n9C999mrr mnd id9ntity by biock numbmr) 

wind Tunnel Wall Corrections 
Vortex Lattice Method 


20. abstract ^Conf/nu« on rovorao «/d« // nmc999mry mnd idmntify by block mjmbmr) 

A computer program was developed to obtain the wind 
tunnel wall corrections for wing angle of attack, induced drag, 
and pitching moment in incompressible flow. The vortex 
lattice method is used for computation of these correction 



DD,j°n“j 1473 edition of I NOV «s is obsolete UNCLASSIFIED 

S/N 0t02-014- 6601 | 



1 



1CCU«ITY CLAlJiriCATlON OF THIS ^AOC -D### Mntmr^d) 



UNCLASSIFIED 



f li C u TV CL ASSI A T»On This g4GCr^#n Orntm 

(20. ABSTRACT Continued) 

factors. The program can be applied to wind tunnels of 
arbitrary cross-sectional shape, and wings of any desired 
planform, subject to the constraint of straight leading 
and trailing edges. 



DD Form 1473 
1 Jan 73 

5/N 0102-014-6601 



UNCLASSIFIED 



^ security cl AlllElCATIOM THIS P kGtfWh*’' Cnt»r»rf) 



liiMf 




Approved for public release; distribution unlimited. 



Wind Tunnel Wall Corrections for Arbitrary 
Planforms and Wind Tunnel Cross-Sections 

by 

Chester Arthur Heard 
Lieutenant, United States Navy 
B.A., Southwestern at Memphis, 1969 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 

from the 

NAVAL POSTGRADUATE SCHOOL 
June 1977 



f^woxuBWRr 

’OSTGMOUATESCHOOl 



ABSTRACT 



A computer program was developed to obtain the wind 
tunnel wall corrections for wing angle of attack, induced 
drag, and pitching moment in incompressible flow. The 
vortex lattice method is used for computation of these 
correction factors. The program can be applied to wind 
tunnels of arbitrary cross-sectional shape, and wings of 
any desired planform, subject to the constraint of straight 
leading and trailing edges. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 12 

II. METHODS USED 15 

A. GENERAL 15 

B. VORTEX LATTICE METHOD 

C. AIRFLOW REPRESENTATION ^0 

D. DETERMINATION OF WALL CORRECTIONS 

III. RESULTS AND DISCUSSION 30 

A. THE COMPUTER PROGRAM 30 

B. GROUND EFFECT 31 

C. WIND TUNNEL APPLICATION 34 

IV. CONCLUSIONS AND RECOMMENDATIONS 43 

APPENDIX A: ASSUMPTIONS 44 

APPENDIX B: DEVELOPMENT OF INDUCED 

VELOCITY EQUATIONS 45 

APPENDIX C; COMPUTER PROGRAM INPUT 50 

COMPUTER PROGRAM 51 

LIST OF REFERENCES 74 

INITIAL DISTRIBUTION LIST 75 



5 



LIST OF FIGURES 



1. Wind Tunnel Geometry 

2. Wing Geometry 

3. Ground Effect Ratio Plotted as a Function 

of Wing Height Above the Ground 

4. Wing Span Load Distribution 

5. Induced Drag Distribution 

6. Angle of Attack Correction Factor Plotted 

as a Function of Span Ratio 

7. Angle of Attack Correction Factor Plotted 

as a Function of Sweep Angle 

8. Velocity Induced by a Vortex Segment 



18 

22 

32 

35 

39 

40 

41 
46 



6 



LIST OF TABLES 



Ground Effect Factor Variation With 
Sweep Angle 

Wing Parameter Variation with Span 

Wing Parameter Variation with 
Sweep Angle 



LIST OF SYMBOLS 



[AERO] 


wing aerodynamic influence coefficient matrix 
at the 0.75 chord position 


AR 


2 

Wing aspect ratio = b /S 


b 


Wing Span 


C 


Wind tunnel cross-sectional area 


c 


Local wing chord, streamwise direction 


^ave 


Average wing chord, (s/b) 




Wing span loading, (H/q) 


^L 


Wing lift coefficient = Lift/qS, dimensionless 


^L /^L 
a 

g Fa 


Ground effect ratio 




Pitching moment coefficient = Moment/qSc 


S 


Root chord 


S 


Tip chord 


s. 

1 


Wing induced drag coefficient 




Spanwise induced drag distribution (d/q) 


d 


Wing induced drag distribution, drag/unit span 


[DAERO] 


Wing aerodynamic influence coefficient matrix 
at the wing quarter-chord 


[GAMAL] 


Wall vortex loop strength 


h 


Wing height above the ground 


H 


Wind tunnel height 


1 


Wing running span load, lift/unit span 


NS 


Number of horseshoe vortices per semispan 


q 


Dynamic pressure 



8 



[ROL] 


Wall aerodynamic influence coefficient matrix 


[ROT] 


Influence coefficient matrix for wall-on-tail 


[ROW] 


Influence coefficient matrix for wall on wing 


S 


Wing area 


TR 


Wing tape ratio, C /C„ 


V 


Freestream velocity 


w 


Induced downwash velocity 


W 


Wind tunnel width 


3c/4 


Induced downwash angle at 0.75 chord position 


[WOL] 


Influence coefficient matrix for wing-on-wall 


[WOT] 


Influence coefficient matrix for wing-on-tail 


X 


Axial coordinate 


Y 


Vertical coordinate 


Z 


Spanwise coordinate 


a 


Angle of attack 


6 


Wind tunnel wall interference factor 


r 


Vortex strength 


£ 


Induced downwash angle at the tail 


n 


2 Z 

Wing spanwise stations (-^) 


*^TL 


Tail efficiency factor 



Subscripts : 



FA 


free air 


g 


ground 


1 


wall vortex loops 


T 


wind tunnel 



9 



TL 

w 



tail 
wing 
9( )/3a 



Matrix Notation: 



[ ] 

{ } 

r' j 

[ ] 



-1 



Rectangular matrix 
Column matrix 
Diagonal matrix 
Inverse of a square matrix 



10 




I 



D 





ACKNOWLEDGMENT 



My sincere gratitude goes to Professor L. V. Schmidt 
for his invaluable assistance throughout the course of this 
work and to my wife Anne for her patience and understanding. 



11 



I. 



INTRODUCTION 



In the process of wind-tunnel testing the experimental 
procedure itself constrains and alters the flow past the 
wing under examination. This modified flow induces errors 
which must be corrected if accurate results are to be 
obtained. The major interference effects in subsonic wind 
tunnels result from the constraining influence of wind 
tunnel walls, disturbances from model supports and measuring 
equipment, and turbulence and non-uniformity in the airflow 
itself . 

Russell [Ref. 10] treats the problems resulting from 
model support interference and introduces a method for 
reduction of this effect. The differential equations 
governing flow within a wind tunnel are the same as those 
in free air. However, the presence of the wind tunnel walls 
changes the outer boundary conditions. The result is a 
series of perturbations of the velocity-potential gradient 
normal to the model surface. The effect is pervasive in 
that almost all of the measured aerodynamic quantities must 
be corrected. 

Generally, Prandtl is acknowledged as having established 
the foundations for research on wind-tunnel wall interfer- 
ence. His lifting-line theory and concept of trailing vor- 
tices opened the doors to much experimental investigation 
in this area. The resulting studies considered two-dimensional 
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and three-dimensional effects, open and closed wind tunnels, 
and various geometrically shaped test sections. 

References 1 and 8 offer a comprehensive treatment of 
wind-tunnel wall interferences. Much tabular and graphical 
data are presented for different tunnel shapes and wing 
spans. Pankhurst and Holder included much data on octagonal 
wind tunnels, while Pope concentrated more on a uniformly 
loaded wing. 

With the advent of sweptback, slender wings designed 
for high-speed flight, the lifting line model was no longer 
suitable. In these situations, the wing lifting system 
may be represented by distributed horseshoe vortex elements 
[Refs. 6 and 7]. In Ref. 9, data, graphs, and formulas are 
easembled from various sources and systematically analyzed 
for various types of wings in wind tunnels . Although the 
mathematical formulations are tedious, the resulting graphs 
and tables summarized very well the data obtained from each 
theoretical treatment. 

Even with all this information available, it is still 
an arduous task to approach tabulated data with a specific 
wind tunnel shape and geometric configuration, and hope for 
a solution to the wall interference problem. It is the 
purpose of this paper to develop a computerized, closed- 
form solution to this problem. The design goal was to 
develop a digital computer code that would be applicable to 
wind tunnels of arbitrary cross-section and wings with 
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varying planforms. Such versatility would enhance the 
accuracy of test data and provide a valuable, accurate 
tool in wind-tunnel testing. 
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II. METHODS USED 



A. GENERAL 

Wind tunnels provide an invaluable tool in the design 
and performance prediction of wings and wings-and-bodies . 
However, the wind-tunnel testing environment introduces 
its own set of problems. The interference action from the 
walls of the wind tunnel alters the development of induced 
velocities by the wing lifting system. As a result, downwash 
formation is suppressed and the wing aerodynamic character- 
istics are those of a wing having greater aspect ratio. The 
accurate determination of wind tunnel wall interference is, 
therefore, an important consideration in wing experiments. 

The classical wall correction factor, 6, for the case 
of a single horseshoe vortex representing the complete wing 
is defined by; 



Aa = 




( 1 ) 



where 



A . -1 w 

Aa = tan ^ 

Using small angle approximations, the tangent is approxi- 
mately equal to the angle. The wall-correction factor may 
be expressed in terms of wing circulation, vortex span and 
average induced downwash velocity by: 



15 



(5 



( 2 ) 



w c 

2 b r 

V 

Calculation of 6 is dependent on the determination of 
the downwash velocity induced by the wind tunnel walls. 
Classical methods involve the utilization of an image 
system of vortices which cancel the effect of the wing 
lifting system at the wall. The boundary condition of no 
flow through the wall is, therefore, preserved. This system 
can provide an accurate, closed-form solution to the ground 
effect problem. However, when four walls are present, each 
wall must be represented by an image system which affects 
other image systems, as well as the model itself. The 
result is a series of infinite systems. Truncation of this 
system introduces inaccuracies of unknown magnitude and 
significance. 

This method of images has been applied to tunnels of 
rectangular cross-section with success. However, it is 
inadequate when arbitrary tunnel cross-sections are 
encountered. Add to this a wing with sweep and taper, and 
the problem becomes significantly altered. 

Therefore, a new method was sought to solve this problem 
specifically for the NPS 3.5 X 5 ft. octagonal wind tunnel. 
The technique chosen was presented by Joppa in Refs. 4 and 
5. This approach lends itself very well to computer solution 
and can be applied to wind tunnels of any cross-sectional 
shape . 
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In this study, an extension of Joppa's vortex lattice 
method was made to allow the treatment of wings with varying 
shape, sweep, and aspect ratio. Such treatment requires a 
modification of the single horseshoe vortex representation 
of the wing lifting system. The method of Weissinger as 
presented in Ref. 6 provides an accurate means of representing 
variable planforms. 

B. VORTEX LATTICE METHOD 

Since the method of images was considered to be limited 
in application, it was necessary to find another, more versa- 
tile method of representing the wind-tunnel walls. The 
problem is essentially one of finding a vorticity distribu- 
tion which will accurately model the wind-tunnel walls. This 
vortex system must satisfy the boundary condition of no 
flow through the walls. 

In Ref. 4, Joppa presents a method which accurately 
models the wind-tunnel walls. This was accomplished by 
representing the tunnel walls with a lattice network of 
rectangular vortex loops. This vortex system lies in the 
plane of the tunnel walls and satisfies Helmholtz's theorem 
that a vortex filament can neither begin nor end in the 
flow. Fig. 1 shows a representation of this system as 
applied to the NFS 3.5 by 5 ft. wind tunnel. 

Each vortex segment in a given rectangle has vortex 
strength F^, and has a corresponding control point located 
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Figure 1. Wind Tunnel Geometry 
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ire. 



in the center of each rectangle where the no-flow boundary 
condition is satisfied. The last vortex lattice rectangles 
located far downstream are represented by horseshoe vortices 
trailing to infinity. The reasoning is that far downstream, 
where the influence upon the wall vortex lattice array of 
the bound vortex elements on the wing is small relative to 
the trailing vortex elements, one finds that the vortex 
strengths of the lattice terms do not change in the stream- 
wise direction. Consequently, by a straightforward cancella- 
tion argument, one can visualize replacing the vortex lattice 
rectangles by horseshoe vortices trailing to infinity. 

The mathematical presentation of the theory is based on 
the Biot-Savart law describing the velocity induced at a 
point by a vortex segment. Joppa in Refs. 4 and 5 follows 
the mathematical formulation beginning with the Biot-Savart 
law. This development is summarized in Appendix B. 

The general matrix representation of this method centers 
around the boundary condition of no flow through the wall. 

As mentioned previously, this condition is satisfied at a 
control point located at the center of each vortex loop. 
Therefore, it follows that the velocity induced at each 
control point by the vortex loops must be exactly cancelled 
by the velocity created by the wing lifting system. This 
is presented below in matrix representation 

[ROL] {r^} = [WOL] {r^} (3) 
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which allows a solution of vortex lattice strengths in 
terms of the wing vortex strengths. 

{r_} = [ROL]"^ [WOL] {r } (4) 

ij W 

Note that the [ROL] matrix, although possibly quite large, 
is square and hence may be inverted. 

[ROL] represents the wall influence coefficient matrix. 
It is obtained by calculating the downwash induced at each 
wall control point by a unit value of each vortex loop. 
Therefore, ROL^^ is the downwash at wall control point "i" 
induced by vortex loop "j". [WOL] is the influence coeffi- 
cient matrix describing the wing effect on each loop con- 
trol point. WOL^j is the downwash produced at control point 
th 

"i" by the j wing horseshoe vortex. 

By inverting [ROL] and premultiplying with [WOL] , the 
strength of the wall vortex loops can be found in terms 
of the wing circulation. The inversion of [ROL] can be 
very time-consuming, since it may be a very large matrix. 
However, by using symmetry, the order of this matrix can 
be reduced by a factor of four. This requires the wing to 
be placed on the horizontal and vertical centerline of the 
wind tunnel. 

C. AIRFOIL REPRESENTATION 

In Ref. 5, Joppa stated that the wing can be repre- 
sented by a lifting system more complicated than a single 
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elementary horseshoe vortex. In order to model swept wings 
of arbitrary planform, a more definitive lifting system 
model was employed. Gray and Schenk, Ref. 6, present a 
modified Weissinger method designed to find the load dis- 
tribution of wings with arbitrary planform. The wing is 
represented by a series of horseshoe vortices. The midspan 
of each horseshoe vortex bound element is located at the 
wing quarter-chord position. The condition of flow tangency 
is satisfied at the wing three-quarter-chord point in the 
center of each horseshoe vortex. Figure 2 illustrates the 
configuration chosen. 

Pearson [Ref. 7] modified this technique by allowing 
the bound vortices to be swept, and including chord-wise 
horseshoe vortices. For this presentation it was felt that 
unswept bound vortices would accurately represent the lifting 
system of a swept wing. In addition it was felt that little 
deterioration of accuracy would result if the chordwise 
system of bound vortices was combined into a single bound 
vortex at the quarter-chord. Figure 2 illustrates the 
configuration chosen. 

An important concept in this representation of the wing 
is that of flow tangency at the three quarter-chord. The 
downwash angle at the wing can, therefore, be determined 
by placing a control point at this position in the center 
of each horseshoe vortex. The rationale behind the three- 
quarter chord flow tangency condition is a product of 
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Figure 2. Wing Geometry 
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( 



I 




two-dimensional flow considerations. Extension of this 
technique to the three-dimensional wing has been shown to 
produce experimentally confirmed results. However, it is 
still unclear why this technique works so well. 

The wing aerodynamic influence coefficient matrix was 
then developed. [AERO] represents the downwash induced at 
each 0.75 chord wing control point by a unit value of wing 
horseshoe vortices . The matrix formulation for induced 
downwash formation in free air is presented below; 

{w} = [AERO] {r } (5) 

w 

Dividing the above expression by freestream velocity, one 
obtains an expression for 0.75 chord control point induced 
downwash angles. The flow tangency requirement equates 
this column matrix to the wing angle of attack column 
vector. 



= [AERO] {^} (6) 

3c/4 



and 



{a} = 



{-} 



3c/4 



I 2r 

It should be noted here that — is equal to from 



the Kutta-Joukowski Law. Substitution into equation 6 
yields the following results: 



23 




3c 

4 



0.5 [AERO] {-} 

^ FA 



( 7 ) 



£ 

Solving equation 7 for {— }, one can obtain the wing span 

% 

load distribution. Since — is equivalent to the sectional 

<3 

I 

value of C^c, integration of the — distribution vs. span 



q 



provides a direct measure of the wing lift coefficient, . 
When one assumes a uniform angle of attack value of one 
radian along the span, the resultant integration of the 
span load solution will yield the wing lift curve slope. 

It should be borne in mind that these techniques are 
predicated upon principles from thin airfoil theory and 
in the limit for an infinite aspect ratio wing will yield 
a wing lift curve slope of 27r(rad ^) . 

D. DETERMINATION OF THE WALL CORRECTIONS 

The classical wall interference factor 6 , which includes 
a ratio of wing to tunnel cross-section area, is generally 
found by representing the wing as a single horseshoe vortex. 
Since the present formulation utilizes ten horseshoe vor- 
tices per semi-span, 6 in a modification can best be found 
by continuing our matrix algebra. 



a 



FA 




C 



L 



a 



FA 



C 



( 8 ) 



L 



a. 



FA 
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a 






By assvuning that the lift coefficient is the same in 
both tunnel and free air, the following equations result: 



+ 6 C, 



(9) 



a 



FA 



ot. 



Cancelling the lift coefficient and solving for 5 



6 = 



( 10 ) 



a 



FA 






The free air lift slope was found in the previous section. 
The wing lift curve slope for the wing in the wind tunnel 
can be found by a simple extension of this same method. 

The total downwash at the wing control points is equal to 
that induced by the wing added to that created by the tunnel 
walls. The effect produced by the tunnel walls is repre- 
sented by the vortex rectangle-on-wing influence coefficient 

th 

matrix [ROW] . ROW^ . is the downwash induced at the i wing 

th 

control point by the j wall vortex loop. Multiplying this 
matrix by the vortex loop strength, one obtains the total 
wall effect. 



{w} = [AERO] {r } + [ROW] {r_} 

W li 



( 11 ) 
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Substitution of equation 4 for {F } forms the following 

Xj 

set of equations: 

{w} = [AERO] {r } + [ROW] [ROL]~^ [WOL] {F } 

w w 



curve slope can be determined. It should be noted that the 
correction factor obtained from this process is different 
from the classical 6, since the correction factor computed 
above already contains the wing area to tunnel area ratio 
influence. Therefore, the correction to angle of attack 
can be obtained by a direct application of this 5. 



In most literature, the correction for induced drag 
is obtained by application of the same 5 as that used for 
angle of attack. 




3c 

4 



0.5 [[AERO] + [ROW] [ROL]"^[WOL] ] {X} (12) 

^ rp 



T 



By establishing one radian angle of attack, solving 
for — , and integrating the result, the wind-tunnel lift 



Aa 



6 C 



L 



(13) 



AC 



5 C 



2 



(14) 



D. 

1 
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However, in this presentation, the induced drag correc- 
tion is calculated in a different manner because of added 
available aerodynamic information. The spanwise induced 
drag distribution in free air is given by the following 
formula : 




['^^] [DAERO] 
^ FA 




(15) 



[DAERO] is the wing aerodynamic influence coefficient 
matrix computed with the wing control point now located 
at the quarter-chord position. The wall effect can similarly 
be found by including additional matrix terms for the wall 
effect on the wing quarter-chord position. 

{-} = [[DAERO] + [ROW] [ROL] [WOL] ] {-} (16) 

q ^ q .p q T 



Integration of equations 15 and 16 result in the free 
air and wind tunnel induced drag, respectively. The induced 
drag correction factor can then be found in the following 
manner . 



C 





C 



D. 

1 



T 
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K 



FA 



C /C 

D. 



K„ 






K 



^FA 



AC^ = K C, 



( 17 ) 



When testing a model in the tail-on configuration, 
the presence of the wind tunnel walls serves to modify the 
downwash angle induced at the tail. This altered flow 
field can have a significant effect on the pitching moment 
coefficient. The following increment in pitching moment 
coefficient is generated by the horizontal tail: 



AC 






(n 



TL 



TL 



) ct 



w 



a 



TL 



TL 



( 18 ) 



The term in parentheses is approximately equal to the 

3C 

stabilizer effectiveness, „ . The tail angle of attack 

O 0 . 

1 

can be estimated as follows: 



a 



TL 



a - £ + 6 . 

1 



( 19 ) 



At a given value of wing lift coefficient an incremental 
change in downwash angle at the tail due to the wall effect 
can be computed : 
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( 20 ) 



T FA 



Therefore, the correction to pitching moment coefficient is 
given by the following relation: 



AC 



m 



= C . Ae 
mi 



( 21 ) 



Since the downwash angle e is equal to ^ at the tail 
position, £ for free air and for the wind tunnel can be 
found from the following relations, respectively: 



'FA 



0,5 [WOT] {^} 



FA 



= 0.5[[WOT] + [ROT] [GAMAL] ] {-} (22) 

T q ^ 

[WOT] is the wing effect on the tail position and 
[ROT] [GAMAL] represents the wall effect on the tail. 
Equation 20 can then be used to find Ae. The pitching 
moment correction, as evaluated in equation 21, is thus 
dependent on being given an experimental value for the 
tail effectiveness parameter. 
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III. RESULTS AND DISCUSSION 



A. COMPUTER PROGRAM 

The computer program developed consists of two major 
subroutines. Subroutine Loop calculates the downwash at 
any point due to a vortex rectangle. Subroutine HSHOE 
computes the downwash at any arbitrary point due to a 
horseshoe vortex. HSHOE is used for both the wing horseshoe 
and those of the wall system located far downstream in the 
tunnel walls. The program offers two basic options. One 
option establishes whether ground effect or wind-tunnel 
wall corrections are to be calculated. The other option 
indicates whether the tail-on or tail-off case is to be 
considered. 

Each computer run offers the inclusion of many wings 
of different span, aspect ratio, taper ratio, and sweep. 

The wing parameters are read as input data after the inver- 
sion of the large wall-on-wall coefficient matrix. There- 
fore, with a given tunnel configuration, a variety of 
different wings can be included with a minimum of computer 
time consumed. 

Integration of the wing span loads in accomplished by 
a subroutine obtained from the NPS source library. It 
should be noted that in the integration scheme the assumption 
is made that the wing loading at r\ = 0.0 is the same as that 
for Ti = .05b. Any inaccuracy resulting from this assumption 
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appears to be negligible. Appendix C describes the necessary 
input to the program. 

B . GROUND EFFECT 

Originally, ground effect calculations were conducted 
merely to check the accuracy of the vortex lattice method. 

As previously stated, the method of images provides a 
closed-form solution to the ground effect problem. There- 
fore, this provided an excellent opportunity for a conver- 
gence check. As an object of academic interest an option 
for ground effect computation was retained in the final 
computer program. 

The results of the convergence check showed that the 
vortex lattice method could very accurately compute the 
velocity induced at a point by the presence of a wall. By 
applying the theory previously developed, the lift curve 
slopes for free air and ground effect were computed. Results 
were then compared with those presented in Etkin [Ref. 3] . 

The free air lift curve slope was found to correlate extremely 
well for both straight and sweptback wings. 

Furthermore, in Ref. 3, Etkin defines a ground effect 
ratio as the ground effect lift curve slope divided by that 
in free air. The results are tabulated in Table 1, and 
Figure 3 compares the results obtained from the computer 
program with those presented in Etkin. As can be seen, 
agreement is very favorable. 
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GROUND EFFECT RATIO 




b 



Figure 3. Ground Effect Ratio Plotted as a Function 
of Wing Height Above the Ground 
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.25 Chord Sweep Angle, deg 
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Ground Effect Factor Variation with Sweep Angle 



Several computer runs were conducted at different wing 
heights above the ground. For each run, sweep was varied 
from 0 to 40 degrees. Results are tabulated in Figure 4. 

It can be noted that the ground effect ratio initially 
increases with sweepback but returns to the uiiswept value 
at about 30 degrees. 

As a note of caution, care should be exercised when 
the wing gets very close to the ground. It is necessary 
to greatly decrease the vortex loop size for accurate 
results. This procedure seems to provide a much more finely 
defined ground representation. Otherwise, calculations 
will result in inordinately high values. 

C. WIND TUNNEL APPLICATION 

A convergence check of the wind tunnel wall correction 

factor was obtained by applying the program to a 3 X 5 ft 

rectangular tunnel. The correction factor obtained was 

g 

converted to the classical 6 by including the factors ^ 
and 57.3 (radians to degrees). The result of 0.1216 was 
compared to that presented by Joppa in Ref. 5. The wing 
geometric configuration presented in Ref. 5 consisted of a 
single horseshoe vortex. The value of 0.127 from Ref. 5 
correlates very well with the results of this study. 

Tables 2 and 3 present the data obtained from compu- 
tations with the NPS 3.5 X 5.0 ft octagonal wind tunnel. 
Computer runs were conducted using wings of varying span 
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Figure 4. Wing Span Load Distribution 
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wing Parameter Variation With Span 





in 


in 


o 


o 


CTi 


VD 








o 


CO 


CO 




fH 


O 








in 


CN 




in 


CTi 


fH 








r^ 


CTi 


r- 


VD 


VD 


O 






o 


• 


• 


• 


• 


• 


• 








m 


m 






in 










fH 


CO 


CO 


in 




VD 








O 


CO 


o 


ro 


00 


O 










m 


CN 


CN 


CTi 


fH 










«H 


00 




VD 


O 






in 


• 


• 


• 


• 


• 


• 






ro 


ro 










VD 










CN 


ID 


cn 


in 


VD 










ID 


o> 


r- 


00 


O 








CTi 


O 






o 


fH 








o 


cn 


CO 






O 






o 


• 


• 


• 


• 


• 


• 




0) 


ro 
















fH 














C3^ 




tT> 






in 


CN 


o 


in 


VD 




C 




in 


in 


in 


cr> 


r- 


O 




C 




o 


cn 


CN 


fH 


o 


fH 








CN 




cri 


00 




O 




a 


in 


• 


• 


• 


• 


• 


• 




0) 


CN 
















0) 














fH 




15 






CN 


in 






r- 




CO 




O 


fH 


CO 


CTi 


o 


o 








<Tk 


m 


in 




tH 


fH 




x: 




CN 


in 


cr> 


00 


r- 


O 






o 


• 


• 


• 


• 


• 


• 




•H 


CN 












CN 








<y> 


ID 


VD 




CTi 


r- 


00 


a 




r- 


ID 


O 


CTi 


CN 


o 




0 






<T» 


00 


VD 


fH 


fH 


0) 


*H 




ro 


in 


(Ti 


CO 


r- 


O 


fH 


4J 


in 


• 


• 


• 


# 


• 


• 


XI 


ai 
















a 


•H 














CN 


Eh 


H 




vn 


r-* 




CTi 


rH 


r- 




(U 




fH 




CN 


CTi 




o 




> 




CO 


cn 


CT> 




fH 


fH 








ro 


KD 


<Ti 


00 




O 




U 


O 


• 


• 


« 


• 


• 


• 




0) 


fH 
















-P 














00 




<D 




cr> 


in 


CO 






r- 




S 




CN 


r- 




CN 




o 




CC3 








CJ> 


00 


fH 


iH 




H 




cn 


VD 


cr> 


CO 


r- 


O 




cd 




• 


• 


• 


• 


• 


• 




cu 


in 
















Cn 




in 




CN 




VO 


r- 




C 






o 






00 


o 




•H 




in 


cr> 


m 


CN 


o 


fH 




:2 




CN 




<Ti 


CO 




O 








# 


• 


• 


• 


• 


• 






o 












C 
















a 


0 
















0 


*H 






a 










•H 


4J 






0) 












O 






0) 










O 


0) 
















<D 


H 






CQ 










U 


H 














U 


0 








< 








0 


U 






J-l 


Cli 


Eh 


< 




U 








rr 0 




S 


IH 


Eh 




Q 






\-C 






Q 


Q 


d 


U 






O U 


U 


u 


U 


U 


<3 


<1 







37 



and sweep angle. Figures 4 and 5 illustrate the span loads 
calculated for a typical wing. It is interesting to note 
that as sweep angle increased, negative induced drag appears 
on the outboard wing stations. Therefore, total induced 
drag is seen to steadily decrease with increasing sweep 
angle. 

As sweep angle is increased, the free air lift curve 
slope decreases. Correspondingly, the corrections to 
angle of attack and induced drag are reduced. Figure 6 
graphically displays the data for a straight wing. In the 
past, the correction factors for the NFS 3.5 X 5.0 ft 
octagonal wind tunnel were obtained by assuming an elliptic 
cross-section. The correction factors presented in Ref. 1 
for elliptic cross-section are 0.610 for angle of attack 
and 0.0106 for induced drag. In comparison, the vortex 
lattice program computes values of 0.7086 and 0.01075, 
respectively. It is interesting to note that the induced 
drag correction is almost identical, whereas the correction 
to angle of attack is slightly higher for the octagonal 
wind tunnel. Figure 7 illustrates the angle of attack 
correction factor with varying tunnel span ratios. The 
correction factor can be seen to steadily increase with the 
span ratio. 

It is an interesting observation that the value of 
induced drag for a wing span of 5.0 ft was almost zero. 

This compares very well with airfoil theory in that a wing 
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Figure 5. Induced Drag Distribution 
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Figure 6. Angle of Attack Correction Factor Plotted 
as a Function of Span Ratio 



40 




SWEEP ANGLE fJ.qj 



Figure 7. Angle of Attack Correction Factor Plotted 

As A Function Of Sweep Angle 
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which spans the wind tunnel represents an approximation to 
the two-dimensional wing. Induced drag for a two-dimensional 
wing should equal zero. In addition, the wind-tunnel lift 
curve slope is within 0.0006 of the theoretical two-dimensional 
value of 2tt rad 

A sample computer run was made for the moment correction 
factor. The tail quarter-chord was placed two feet behind 
the wing bound vortex, on the tunnel centerline. The 
correction factor obtained for a straight wing of aspect 
ratio 6.0 was -0.2297. For a wing with sweep angle of 
30 degrees and taper of 0.3, a value of -0.1887 was calcu- 
lated. The above correction factors must be multiplied 
by the stabilizer effectiveness parameter, using consistent 
angular measures, in order to find the pitching moment 
correction. 
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IV. CONCLUSIONS AND RECOriMENDATIONS 



It is felt that the vortex lattice computer program 
as presented in this paper offers a very reliable and 
flexible tool for the determination of wind-tunnel wall 
corrections. It is capable of treating wings of any span, 
taper ratio, aspect ratio, and quarter-chord sweep, subject 
to the constraint of straight leading and trailing edges. 

The geometric shape of the wind tunnel cross-section can 
similarly be varied at the discretion of the program user. 
Resultant accuracy appears to be excellent. However, cau- 
tion should be exercised in determining the vortex loop size. 
As the wing is positioned closer to the walls, a smaller, 
more definitive vortex loop representation must be used. 

A logical extension of this program would be the inclu- 
sion of wings with an antisymmetric load distribution. 
Furthermore, slight modifications could incorporate the 
ability to consider the wall effects on other stability 
derivatives. In its present form, this program offers 
the wind tunnel user a versatile tool for determination of 
wind tunnel wall corrections. 
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APPENDIX A 



The general assximptions used in the compilation of 
this program are listed below. 

1. The two-dimensional lift-curve slope is 2ir rad 

2. The wing load distribution is symmetrical. 

3. Aeroelastic effects were not considered. 

4. The wing is located on the horizontal and vertical 
centerline of the wind tunnel. 

5. The leading and trailing edges of the wing are 
straight. 
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APPENDIX B 



DEVELOPMENT OF INDUCED VELOCITY EQUATIONS 

The equations required to solve for induced velocities 
consist of two main types. One equation solves for the 
velocity induced by an arbitrarily oriented vortex segment, 
whereas the other basic equation applies to horseshoe vortex 
elements. These equations have been incorporated in sub- 
routines LOOP and HSHOE, respectively. The development of 
these equations presented below follows closely that found 
in Ref. 4 and 5. 

According to the Biot-Savart law, the velocity induced 
at a point p is found from: 



w 



-7— r- (cos B, + cos B_) V 
4irh 1 2 



(IB) 



For vortex segments oriented in the Y-Z plane. 



cos B^ + cos B^ 



^1 ^2 2 2 



(2B) 



It is now necessary to consider v, the unit vector 
which establishes direction. From Fig. 8, it can be seen 
that 



_ X S 

V = — — — 

|r^ X s 



(3B) 
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Figure 8. Velocity Induced By A Vortex Segment 
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Noting from Fig. 8 that 



hs = |r^ X S 



(4B) 



then 



hs = + (R^^S^ - R^^S^)2 

+ (R^ S - R^ S 
lx y ly X 



Calculating the induced velocity in the direction of 
outward normal of point P per unit F, 



- V — 

P - p • n 



(5B) 



where the unit normal vector is defined by: 



n = oi + (sin 9p) j + (cos 0p)k 



Therefore, one finally obtains: 



^1 ^2 2 2 
i ^ [S^ - (R, - R-)^] 

8TTR^R2(hs)^ ^ 

(6B) 

' [(sin 9P) + (COS0P) (R^^S^- R^yS^)] 
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■A 




A conditional statement must be included in the program 
to set the downwash equal to zero if hs is equal to zero. 
This will preclude having a zero term in the denominator. 

The trailing elements of a vortex loop can be treated 
with more simplicity. For these segments. 



cos 



cos 




- 



- 



(7B) 



However, it can be seen from Fig. 8 that 



R 



h 



1 

2 



2 



2 2 2 
(xp - xA) + (yA-yp) + (zA-zp) 

2 2 
(yA - yp) + (zA - zp) 



(8B) 



Combining equations 7 and 8 , 



cos Bj^ 




cos B 2 




(9B) 



As a result, the equation for the trailing elements 
of a vortex loop is simply 
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w 



cos B 



1 



+ cos B 



2 



r 




[ (sin 0p) (zA - zp) 



(lOB) 



+ (cos 6p) (yp - yA) ] 



For the horseshoe vortex, the equation for the bound 
segment is identical to equation 6. Similarly, the 
trailing elements can be treated as in equation 11 with a 
slight modification. In this case, the cosine of one 
enclosed angle is equal to 1, since the trailing element 
extends downstream to infinity. The induced velocity 
equation becomes 



1 + cos B 



— [(zA - zp) j + (yp - yA)k] 



w 



r 



4irh 



w 



1 + cos B 



— [ (sin 6p) (zA - zp) 



r 



4irh 



(IIB) 



+ (cos 0p) (yp - yA) ] 
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APPENDIX C 



COMPUTER PROGRAM INPUT 



Card Number 


Column Number 


Description 


Format 


1 


1-72 


Program title 


18A4 


2 


1-5 


MM=Number of Spanwise 
Stations 


15 




6-10 


NN=Number of Downstream 
Stations 


15 




11-20 


Deltax=Loop length 


FIO.O 


3 


1 


NOPT=l for Ground 
Effect 


11 




2 


NTAIL=1 for tail-on 


11 


4 


1-72 


Tunnel y-coordinates 
(one quadrant only) 


8F10.0 


5 


1-72 


Tunnel z-coordinates 
(one quadrant only) 


8F10.0 


6 


1-10 


Span 


FIO.O 




11-20 


Aspect Ratio 


FIO.O 




21-30 


Taper Ratio 


FIO.O 




31-40 


c/4 Sweep (degrees) 


FIO.O 


7 


1-10 


x-Coordinate of tail 
quarter-chord 


FIO.O 




11-20 


z-Coordinate of tail 
quarter-chord 


FIO.O 



Input cards 5 and 6 can be repeated for different 
planforms and tail locations. 
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ELU />NC SLVB TC SOLVE :( L/Q )♦( /JERC ) = ( ALPHA ) 
CALL ELL (AERO, NS, 10) 

CALL SLVB (AER0,NS,10,CLC) 

V»PITE(6,2003) 

kFIlE<6t200A) 
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ELU 4ND SLVB TO SOLVE THE MATRIX EQUATION 
(4LPhA)= (AERC-W)«( L/Q) . 
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ESTABLISH THE AERODYNAMIC INFLUENCE COEFFICIENT MATRIX 
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